Fairness in Academic Course Timetabling' 

Moritz Muhlenthaler Rolf Wanka 

Department of Computer Science 
University of Erlangen-Nuremberg, Germany 
{moritz . muehlenthaler , rwanka}(§cs .fau.de 



Abstract 

We consider the problem of creating fair course timetables in the setting of a university. 
Our motivation is to improve the overall satisfaction of individuals concerned (students, teach- 
ers, etc.) by providing a fair timetable to them. The central idea is that undesirable arrange- 
ments in the course timetable, i.e., violations of soft constraints, should be distributed in a 
fair way among the individuals. We propose two formulations for the fair course timetabling 
problem that are based on max-min fairness and Jain's fairness index, respectively. Further- 
more, we present and experimentally evaluate an optimization algorithm based on simulated 
annealing for solving max-min fair course timetabling problems. The new contribution is con- 
cerned with measuring the energy difference between two timetables, i.e., how much worse 
a timetable is compared to another timetable with respect to max-min fairness. We introduce 
three different energy difference measures and evaluate their impact on the overall algorithm 
performance. The second proposed problem formulation focuses on the tradeoff between fair- 
ness and the total amount of soft constraint violations. Our experimental evaluation shows that 
the known best solutions to the ITC2007 curriculum-based course timetabling instances are 
quite fair with respect to Jain's fairness index. However, the experiments also show that the 
fairness can be improved further for only a rather small increase in the total amount of soft 
constraint violations. 



1 Introduction 

The university course timetabling problem (UCTP) captures the task of assigning courses to a 
limited number of resources (rooms and timeslots) in the setting of a unversity. In this work, we 
consider the problem of creating fair course timetables in the context of a particular variant of the 

*Research funded in parts by the School of Engineering of the University of Erlangen-Nuremberg. A preliminary 
version appeared in: Proc. 9th Int. Conf. on the Practice and Theory of Automated Timetabling (PATAT), 2012, pp. 
114-130. 
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UCTP, the curriculum-based course timetabling (CB-CTT) problem proposed in [IDGMS07II . The 
CB-CTT formulation features various hard and soft constraints which model typical real- world re- 
quirements. The hard constraints are basic requirements, for example, no two lectures may be held 
in the same room at the same time. The feasible solutions to a CB-CTT instance are the timetables 
satisfying all hard constraints. Soft constraints characterize properties of a course schedule which 
are undesirable for the stakeholders. The quality of a feasible timetable is determined by the extent 
to which the soft constraints are violated. A soft constraint violation results in a penalty, and the 
task is to find a feasible timetable such that the total penalty is minimal. Situations may arise how- 
ever, in which a large part of penalty hits only a small group of the stakeholders, who would thus 
receive a poor timetable in comparison to others. In other words, a timetable may be unfair due 
to an unequal distribution of penalty. In this work, we apply faimess criteria to the CB-CTT prob- 
lem. Our goal is to achieve a balance of interests between the stakeholders by imposing faimess 
conditions on the distribution of penalty and thus, to improve the overall stakeholder satisfaction. 

In general, faimess is of interest whenever scarce resources are allocated to stakeholders with 
demands. In economics for example, the distribution of wealth and income and how to measure 
inequality of resource distributions is of major concem, see for example [IFS06II and [ISF97I . In 
computer science, faimess is a central theme in design and analysis of communication protocols 
(see for instance [|BFCYZ02[ IBG92[ IBFTlll IJCH841 IKMT98L IKRTOli IOW041 [SB081 ). In opera- 
tions research, faimess criteria have been applied for example to the aircraft landing problem by 
HSKOSn . In the literature, there is a wealth of different definitions of how to determine the faimess 
of a given resource distribution. For instance, we may consider the total amount of allocated re- 
sources, the outcome for the worst-off stakeholder, the deviation from the mean allocation, and so 
forth. We propose two fair variants of the CB-CTT problem, which differ with respect to the under- 
lying notion of faimess. Our first problem forumation, MMF- CB-CTT, is based on lexicographic 
max-min faimess, to which we will refer to as max-min faimess for brevity. This faimess notion 
often appears in the context of network bandwidth allocation (see for example I1BG92[ ISBOSID . 
Max-min faimess is a purely qualitative measure of faimess, i. e., given two allocations, max-min 
faimess tells us which of the two is better, but not by how much. Our second problem formula- 
tion, JFI-CB-CTT, is based on Jain's faimess index proposed by [,JCH84 1. This faimess measure 
is used in the famous AIMD algorithm by [CJ89| used in TCP Congestion Avoidance. In contrast 
to max-min faimess, it conveniently represents the inequality of a resource allocation as a number 
between zero and one. 

In order to solve the MMF-CB-CTT problem we propose MaxMinFair_SA, an optimization 
algorithm based on simulated annealing (SA). Due to the mild requirements of SA on the problem 
stmcture, the proposed algorithm can easily be tailored to other max-min fair optimisation prob- 
lems. A delicate part of the algorithm is the energy difference function, which quantifies how much 
worse one solution is compared to another solution - a piece of information we do not get directly 
from max-min faimess. We propose three different energy difference functions and evaluate their 
impact on the performance of MaxMinFair_SA on the 21 standard instances from track three 
of the ITC2007 competition ( [IDGS12II '). Our experiments indicate, that the known best solutions 
with respect to the CB-CTT model are quite fair in the max-min sense, but further improvements 
are possible for 15 out of 21 instances and often a considerable improvement of the worst-off 
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stakeholder is achieved. 

The fairness conditions imposed by max-min fairness are rather strict in the sense that there 
is no tradeoff between fairness and the total amount of penalty. In practice however, it may be 
desirable to pick a timetable from a number of solutions with varying tradeoffs between fairness 
and total penalty. Our proposed JFI-CB-CTT problem is a bi-criteria optimisation problem, which 
offers this option. We investigate the tradeoffs between fairness and total penalty for the six stan- 
dard instances from [IDGS12I whose known best solutions have the highest total penalty compared 
to the other instances. Our motivation for this choice of instances is simply that if the total penalty 
of a timetable is very small, then there is not much gain for anyone in distributing the penalty in 
a fair way. Our conclusion regarding this approach is that, although the known best solutions for 
the six instances are already quite fair, we can improve the faimess further at the cost of only a 
small increase in total penalty. For a theoretical treatment of the price of faimess on so-called 
convex utility sets with respect to proportional faimess and max-min faimess, see the recent work 
by hBFT lll. 

The remainder of this work is organized as follows. In Section [2l we will provide a brief 
review of the curriculum-based course timetabling (CB-CTT) problem model as well as max-min 
faimess and Jain's faimess index. In Section |3] we will propose MMF-CB-CTT and JFI-CB-CTT, 
two fair variants of the CB-CTT problem formulation, and in Section HI we will introduce the SA- 
based optimisation algorithm MaxMinFair_SA for solving max-min fair allocation problems. 
Section \5\ is dedicated to our experimental evaluation of the faimess of the known best solutions 
to 21 standard instances from the website by [IDGS121 with respect to max-min faimess and Jain's 
faimess index, and the performance of the MaxMinFair_SA algorithm. 

2 Preliminaries 

In this section, we provide a brief review of the curriculum-based course timetabling problem 
formulation as well as relevant definitions conceming max-min faimess and Jain's faimess index. 

2.1 Curriculum-based Course Timetabling Problems 

Curriculum-based Course Timetabling (CB-CTT) is a particular variant of the UCTP. It has been 
proposed in the course of the second intemational timetabling competition in 2007 (see LMSP^IOJ ). 
and has since then emerged as one of the de-facto standard problem formulations in the timetabling 
community. The central entities in the CB-CTT formulation are the curricula, which are sets of 
lectures that must not be taught simultaneously. Both problem formulations proposed in the next 
section are based on the CB-CTT model. 

CB-CTT problems are NP-hard and a lot of effort has been devoted to the development of 
exact and heuristic methods which provide high quality solutions within reasonable time. A wide 
range of techniques has been employed for solving CB-CTT instances including but not limited 
to approaches based on Max-SAT ( HAANIOII ). mathematical programming ( ULLlOi IBMPRUll ). 
local search (rD S06[|LH1 01), evolutionary computation ( [IABM07II ) as well as hybrid approaches 
( LBDS12.,Mul09ll ). There has been a lot of progress in terms of the achieved solution quality in the 
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recent years. Interestingly, there seems to be no approach which is superior to the other approaches 
on most ITC2007 instances (see the website nDGS121 for current results). 

A CB-CTT instance consists of a set of rooms, a set of courses, a set of curricula, a set of 
teachers and a set of days. Each day is divided into a fixed number of timeslots, a pair composed 
of a day and a timeslot is called a period. A period in conjunction with a room is called a resource. 
Each course consists of a number of lectures, i. e. a number events to be scheduled, is taught by 
some teacher and has a number of students attending it. Each curriculum is a set of courses. For 
each room, we are provided with the maximum number of students it can accommodate and for 
each course we are given a list of periods in which it cannot be taught. A solution to a CB-CTT 
instance is a timetable, i. e. an assignment of courses to resources. The quality of a timetable is 
determined according to four hard and four soft constraints (see HDGMSOVID . 

The hard constraints are the following: 

HI All lectures need to be scheduled and no two lectures of the same course may be assigned to 
the same period. 

H2 No two lectures may be assigned to the same resource. 

H3 Two courses in the same curriculum or taught by the same teacher must be assigned to 
different periods. 

H4 A lecture can only be scheduled in a period that is not marked unavailable for the corre- 
sponding course. 

A timetable that satisfies all hard constraints is called feasible. 

The CB-CTT formulation features the following soft constraints: 

51 RoomCapacity: Each lecture should be assigned to a room of sufficient size. 

52 MinWorkingDays: The lectures of each course should be distributed over a certain minimum 
number of days. 

53 IsolatedLectures: For each curriculum, all lectures associated to the curriculum should be 
scheduled in adjacent timeslots. 

54 RoomStability: The lectures of each course should be assigned to the same room. 

Each violation of one of the soft constraints results in a "penalty" for the timetable. The CT- 
CTT objective function aggregates individual penalties by taking their weighted sum. Detailed 
descriptions of how hard and soft constraints are evaluated and how much penalty is applied for 
a particular soft constraint violation can be found in the report by [iDGMS07 1. Given a CB-CTT 
instance /, the task is to find a feasible timetable such the aggregated penalty is minimal. 
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2.2 Fairness in Resource Allocation 



Fairness issues typically arise when scarce resources are allocated to a number of stakeholders with 
demands. Fair resource allocation has received much attention in economic theory (see for example 
IIFS06II ). but also occurs in a wide range of applications in computer science including bandwidth 
allocation in networks ([BG921) and task scheduling ( URDKOSH ). In many optimization problems 
related to resource allocation, the goal is to maximize the total amount of resources allocated to the 
stakeholders. Fairness in this context means that the distribution of resources over the stakeholders 
is important and that certain resource distributions are preferred over others. 

Consider a resource allocation problem with n stakeholders resources. Each resource allocation 
(admissible solution) induces an allocation vector X — (xi, . . . where each item x,, I < i < n, 
corresponds to the amount of resources allocated to stakeholder i. Typically, a preference for 
certain resource distributions is implicitly or explicitly contained in the objective function. For 
example, the task can be to find allocations maximizing the sum of the individual allocations, the 
mean allocation, the root mean square (RMS), the smallest allocation, and so forth (see [ OgrlO[ 
ISK08II ). When the fairness of an allocation is important we may be interested in improving the 
outcome for the worst-off stakeholders or generally try to allocate resources equally among the 
stakeholders. Max-min fairness is a notion of fairness that favours better outcomes for the worst- 
off stakeholders. It has received attention in the area of network engineering, in particular in the 
context of flow control rBF CYZ02[ IKRTOll [SB081 IZLCJ12II . Various inequality measures have 
been proposed such as the Gini index proposed by [|Gin2 n and Jain's fairness index proposed by 
ILJCH84']. Generally, a highly unequal distribution of resources is considered unfair. Our evaluation 
of faimess in academic course timetabling focuses on max-min faimess and Jain's fairness index. 

Our evaluation of faimess in academic course timetabling focuses on the two faimess criteria 
max-min faimess and Jain's faimess index. 



Max-min Fairness. Max-min faimess can be stated as iterated application of Rawls's Second 
Principle of Justice by [|Raw99ll : 

"Social and economic inequalities are to be arranged so that they are to be of greatest 
benefit to the least-advantaged members of society." (the Difference Principle) 

Once the status of the least-advantaged members has been determined according to the difference 
principle, it can be applied again to everyone except the least-advantaged group in order to maxi- 
mize the utility (in the economic sense) for the second least-advantaged members, and so on. The 
resulting utility assignment is called max-min fair. A max-min fair utility assignment implies that 
each stakeholder can maximize his/her utility as long as it is not at the expense of another stake- 
holder who is worse off. Thus, a max-min fair allocation enforces an efficient resource usage to 
some extent. A max-min fair resource allocation is Pareto-optimal. 

In order to define max-min faimess more formally, we introduce some notation. Let X be an 
allocation vector. We generally assume that each entrie of X is a nonnegative real number. By X 
we denote the vector containing the entries of X arranged in nondecreasing order. Similarly, let 
X be a vector containing the entries of X in nonincreasing order. For allocation vectors X and Y 
we write X <mm 1^ if X is at least as good as Y in the max-min sense. For maximization problems 
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such as bandwidth allocation this is the case iff Y ^ig^ X, where :<ig-^ is the usual lexicographic 
comparison. For minimization problems such as the fair course timetabling problems proposed in 
Section [3l X :<mm iff ^ :^iex Y. Let 5 be a solution to an instance / of an optimization problem 
and X be the allocation vector induced by s. Then s is called max-min fair, if for any other solution 
s' to / we have X <mm Y , where Y is the allocation vector induced by s' . Since the allocations 
are sorted, max-min fairness does not discriminate between stakeholders, but only between the 
amounts of resources assigned to them. 

A weaker version of max-min fairness results if Rawl's Second Principle of Justice is not 
applied iteratively, but just once. This means that we are concemed with chosing the best possible 
outcome for the worst-off stakeholder. In the literature, related optimization problems are referred 
to as bottleneck optimization problems ( flEF70[|PZllll ). Note, that in contrast to max-min fairness, 
an optimal solution to a bottleneck optimisation problem is not necessarily Pareto-optimal. In the 
context of practical academic timetabling, the use of bottleneck optimization is hard to justify: 
each stakeholder is guaranteed to be at least as well off as the worst-off stakeholder, but no further 
improvement is considered. 

Jain's Fairness Index. Jain's fairness index is an inequality measure proposed by Jain in [IJCH84II . 
It is the crucial fairness measure that is used in the famous AIMD algorithm by nCJ89ll used in TCP 
Congestion Avoidance. The fairness index J{X) of an allocation vector X is defined as follows: 

( E 4 

J{X) = '-'^ ^ . (1) 

l<i<n 

It has several useful properties like population size independence, scale and metric independence, 
it is bounded between and 1, and it has an intuitive interpretation. In particular J{X) = 1 means 
that X is a completely fair allocation, i. e., the allocation is fair for every stakeholder, and if J{X) = 
l/n then all resources are occupied by a single stakeholder. Furthermore, if J{X) = x% then the 
allocation X is fair for x percent of the stakeholders. 



3 Fairness in Academic Course Timetabling 

Course timetabling problems fit quite well in the framework of fair resource allocation problems 
described in the previous section: A timetable is an allocation of resources (rooms, timeslots) to 
lectures. In this section, we will define two fair versions of the CB-CTT problem formulation 



proposed by f,MSP"'"10l . The first one, MMF-CB-CTT, is based on max-min fairness. Since max- 
min fairness enforces fairness as well as efficiency (maximum utility) at least to some extent, it is 
not a suitable concept for exploring the tradeoff between fairness and efficiency. Thus, we propose 
a second fair variant of CB-CTT called JFI-CB-CTT that is based on Jain's fairness index. 

In order to use the fairness measures mentioned in the previous section, we need determine 
an allocation vector from a timetable. The central entities in the CB-CTT problem formulation 
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are the curricula. Therefore, in this work, we are interested in a fair distribution of penalty over 
the curricula. Depending on the application, a different set of stakeholders can be picked, but 
conceptually this does not change much. Let / be a CB-CTT instance with curricula ci , C2, . . . , 
and let fc be the CB-CTT objective function proposed by P DGMS07 1. which evaluates (S1)-(S4) 
restricted to curriculum c. This means fc determines soft constraint violations only for the courses 
in curriculum c. For a timetable T the corresponding allocation vector is given by the allocation 
function 

A(t) = (/,,(t),/,,(t),...,/,,(t)) . (2) 

Definition 1 (MMF-CB-CTT). Given a CB-CTT instance I, the task is to find a feasible timetable 
T such that A{t) is max-min fair. 

If a feasible timetable has max-min fair allocation vector, then any curriculum c could receive 
less penalty only at the expense of other curricula which receive more penalty than c. This is due 
to the Pareto-optimality of a max-min fair allocation. 

In order to explore the tradeoff between efficient and fair resource allocation in the context of 
the CB-CTT model, we propose another fair variant called JFI-CB-CTT that is based on Jain's 
faimess index proposed by [|JCH84i In order to get meaningful results from the fairness index 
however, we need a different allocation function. Consider an allocation X that allocates all penalty 
to a single curriculum while the remaining k — I curricula receive no penalty. Then J{X) = l/k, 
which means that only one curriculum is happy with the allocation (see [ JCH84II ). In our situation 
however, the opposite is the case: k~l curricula are happy since they receive no penalty at all. The 
following allocation function shifts the penalty values such that the corresponding faimess index 
in the situation described above becomes {k— l)/k, which is in much better agreement with our 
intuition: 

^ ( ~ (/max ~ fci ( '^) 5 fmax ^ fc2 ( 5 • • • 5 fmax ~ fcf. ( ) ? (3) 

with 

fmax = max {/c,(t)} . 

l<i<k 

Definition 2 (JFI-CB-CTT). Given a CB-CTT instance I, the task is to find the set of feasible 
solutions which are Par eto -optimal with respect to the two objectives of the objective function 

F{x) = {f{T),l-J{A'{T))) , (4) 

where f is the CB-CTT objective function from l[DGMS()7^ and J is defined in Eq. ([T]). 

By a similar procedure, other classes of timetabling problem such as post-enrollment course 
timetabling, exam timetabling and nurse rostering can be tumed into fair optimization problems. 
For example, for post-enrollment course timetabling, the central entities of interest are the individ- 
ual students. Therefore, the goal were to achieve a fair distribution of penalty over all students. 
Once an appropriate allocation function has been defined, we immediately get the corresponding 
fair optimization problems. 

Our proposed problem formulations are concerned with balancing the interests between stake- 
holders, who are in our case the students. In practice however, there are often several groups of 



7 



stakeholders with possibly conflicting interest, for example students, lecturers and administration, 
possibilities for extending the problem formulations above to include multiple stakeholders. For 
example, a multi-objective optimization approach may be considered, where each objective cap- 
tures the fairness with respect to a particular stakeholder. When using inequality measures like 
Jain's faimess index for different groups of stakeholders, the inequality values can be aggregated, 
for instance using a weighted-sum or ordered weighted averaging [ Yag88 ] approach. Furthermore, 
max-min faimess or a suitable inequality measure can be applied to the different objectives to 
balance the interests of the different groups of stakeholders. 

4 Simulated Annealing for Max-Min Fair Course Timetabling 

Simulated Annealing (SA) is a popular local search method proposed by I1KGV83 II. which works 
surprisingly well on many problem domains. SA has been applied successfully to timetabling prob- 
lems by [|Kos04|| and IITD96II . Some of the currently known best solutions to CB-CTT instances 
from the ITC2007 competition were discovered by simulated annealing-based methods according 
to the website by [|DGS12 I . Our S A for max-min fair optimization problems shown in Algorithm[r| 
below (algorithm MaxMinFair_SA) is conceptually very similar to the original algorithm. The 
SA algorithm generates a new candidate solution according to some neighborhood exploration 
method, and replaces the current solution with a certain probability depending on i) the quality 
difference between the two solutions and ii) the current temperature. Since max-min faimess only 
tells us which of two given solutions is better, but not how much better, the main challenge in tai- 
loring SA to max-min fair optimization problems is to find a suitable energy difference function, 
which quantifies the difference in quality between two candidate solutions. In the following, we 
propose three different energy difference measures for max-min fair optimization and provide de- 
tails on the acceptance criterion, the cooling schedule, and the neighborhood exploration method 
chosen for the experimental evaluation of MaxMinFair_SA in the next section. 



Algorithm 1: MaxMinFair_SA 
input : Scur: feasible timetable, t^max: initial temperature, t^min: final temperature, timeout 
output: Sbest: Best feasible timetable found so far 

Sbest Scur 
^ ^ i^max 

while timeout not hit do 
Snext ^ neighbor ( Scur) 

'^iPaccept > random then Scur ^ Snext 
if A(Scur) ^mm^(Sbest) then Sbest Scur 

■& ^ next_temperature(t9-) 

end 

return Sbest 
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Acceptance Criterion. Similar to the original S A algorithm proposed by [1KGV831 , algorithm M axMin- 
Fair_SA accepts an improved or equally good solution Snext with probability 1. If Snext is worse 
than Scur then the acceptance probability depends on the current temperature level 'd- and the energy 
difference AE. The energy difference measures the difference in quality of the allocation induced 
by Snext compared to the allocation induced by the current solution Scur- The acceptance probability 
Paccept is defined as: 



where X = A(scur) and Y = A(snext)- In order to fit max-min fairness into the SA algorithm, 
we propose three energy difference measures: AEig^, AEcw, and AEp^. AEjex derives the energy 
difference from a lexicographic comparison, AEcw from the component-wise ratios of the sorted 
allocation vectors and AEp^ from the ratios of the partial sums of the sorted allocation vectors. Our 
experiments presented in the next section indicate that the choice of the energy difference measure 
has a clear impact on the performane of M AxMinFair_SA and is thus a critical design choice. 

For an allocation vector X, let X, denote the ith entry after sorting the entries of X in nonin- 
creasing order. The energy difference AEig^ of two allocation vectors X and Y of length n is defined 
as follows: 



AEiex determines the energy difference between X and Y from the smallest index that determines 
X :<mm Y- Thus, sorted allocation vectors which differ at the most significant indices have a higher 
energy difference than those differing at later indices. In particular, AEig^ has evaluates to 1 if the 
minimum is 1, and it evaluates to 1/n if the minimum is n. 

The energy difference measure AEigj^ considers the earliest index at which two sorted allocation 
vectors differ but not how much the entries differ. We additionally propose the two energy differ- 
ence measures AEcs and AEps which take this information into account. These energy difference 
measures are derived from the definitions of approximation ratios for max-min fair allocation prob- 
lems given by [iKRTOlil . An approximation ratio is a measure for how much worse the quality of a 
solution is relative to a possibly unknown optimal solution. In our case, we are interested in how 
much worse one given allocation is relative to another given allocation. We need to introduce some 
modifications of the definitions by UKRTOin since we are dealing with a minimization problem. 

Note that due to ©, an allocation vector does not contain any positive entries. Let Hxj be the 
smallest value of the two allocation vectors X and Y offset by a parameter 5 > 0, i. e.. 



The offset d is introduced in order to avoid divisions by zero when taking ratios of penalty values. 

The component- wise energy difference AEcw of allocation vectors X and Y is defined as fol- 
lows: 





(5) 



= max{Xi,yi} + 5 



(6) 




(7) 
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Since all entries are subtracted from /ix,y, the ratios of the most significant entries with respect to 
:<mm tend to dominate the value of AEcw Consider for example the situation that Y is much less 
fair than X, say, max{Xi , Yi } occurs more often in X than in Y. Then for a small offset 8 the energy 
difference AEcw{X, Y) becomes large. On the other hand, if X is nearly as fair as Y then the ratios 
are all close to one and thus AEcwiX.Y) is close to zero. 

The third proposed energy difference measure AEp^ is based on the ratios of the partial sums 
a/(X) of the sorted allocation vectors. 

Gi{x) = £ X, . 

The intention of using partial sums of the sorted allocations is to give the stakeholders receiving 
the most penalty more influence on the resulting energy difference compared to AEcw The energy 
difference AEps is defined as: 

A£,.(X,l')= maxfi*£^}-l . (8) 

l<i<n [l-llx,Y-Oi{X) J 

Cooling Schedule. In algorithm MaxMinFair_SA, the function next_temperature updates 
the current temperature level t> according to the cooling schedule. We use a standard geometric 
cooling schedule 

= CC ■ '&max : 

where a is the cooling rate and t is the elapsed time. Geometric cooling schedules decrease the 
temperature level exponentially over time. It is a popular class of cooling schedules which is widely 
used in practice and works well in many problem domains including timetabling problems [1LA87I 
|KAJ94irTD98ll . Geometric cooling was chosen due to its simplicity, since the main focus of our 
evaluation in Section [5] is the performance impact of the different energy difference functions. We 
have made a slight adjustment to the specification of the geometric cooling schedule in order to 
make the behavior more consistent for different timeouts. Instead of specifying the cooling rate a, 
we determine a from -d-max, the desired minimum temperature i9;„m and the timeout according to: 



(■jQ. \ timeout 
'-'max / 

Thus, at the beginning (f = 0) the temperature level is '&max and when the timeout is reached 
{t = timeout), the temperature level becomes d;„in- We chose to set a timeout rather than a maxi- 
mum number of iterations since this setting is compliant with the ITC2007 competition conditions, 
which are a widely accepted standard for comparing results. 

Neighborhood. In our max-min fair SA implementation, the function neighbor picks at random 
a neighbor in the Kempe-neighborhood of Scur- The Kempe-neighborhood is the set of all timeta- 
bles which can be reached by performing a single Kempe-move such that the number of lectures per 
period do not exceed the number of available rooms. The Kempe-move is a well-known and widely 
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Table 1 : Fairness of the known best timetables from [1DGS12II for the ITC2007 CB-CTT instances. 



Instance 


Curricula 


f{sbest) 


J{A'{sbest)) 




compOl 


14 


5 


0.8571 


52,012 




comp02 


70 


24 


0.9515 


4,210,059 




comp03 


68 


66 


0.9114 


l3,lo^9,72,6^5l^4,2^o3^ 

7 3 3 7 7 3 3 7 




comp04 


57 


35 


0.8964 


7,6^,54,42,2,046 




compOB 


139 


291 


0.8277 


4l2,36^35^32^31^30'', 28, 27^262,2514, 


..,2,03 


comp06 


70 


27 


0.9657 


12,72,54,2^,06° 




compOT 


77 


6 


0.9870 


6,0^6 




compOS 


61 


37 


0.9020 


7,6^,54,42,22,049 




comp09 


75 


96 


0.8047 


105,9,710, 66,51°, 4,2,041 




complO 


67 


4 


0.9701 


22,065 




compll 


13 





— 


013 




compl2 


150 


300 


0.9128 


45, 3014 , 28 , 272 , 265 , 25 19 , 224 , 2 16 , 20^ 19, . . 


.,22,03 


complS 


66 


59 


0.8830 


8,7,65,5^42,2^04' 




compl4 


60 


51 


0.9023 


84,7,52,26,04^ 




complB 


68 


66 


0.8495 


103,93,7,64,513,4,2^0^6 




compl6 


71 


18 


0.9176 


72,5^4,06! 




complT 


70 


56 


0.9248 


102,63,59,24,052 




complS 


52 


62 


0.9009 


17,15,14,13,11,10,92,519,22,023 




compl9 


66 


57 


0.9612 


13,7,64,52,4,2^05" 




comp20 


78 


4 


0.9744 


22,0^6 




comp21 


78 


76 


0.8838 


12, 11, 104,9,74,64,512,4,23, l2,045 





used operation for swapping events in a timetable [|BEM+ lOi ILHTOI IMBHSOl ITD981 [TEMOVH . 
A prominent feature of the Kempe-neighborhood is that it contains only moves that preserve the 
feasibility of a timetable. Since the algorithm MaxMinFair_S A only uses moves from this neigh- 
borhood the output is guaranteed to be a feasible timetable. In the future, more advanced neighbor- 
hood exploration methods similar to the approaches in [|DS06[ ILHIOH could be used, which may 
well lead to an improved overall performance of M AXMinFair_S A. 



5 Evaluation 

In this section, we will first address the question how fair or unfair the known best timetables for 
the ITC2007 CB-CTT instances are with respect to Jain's fairness index and max-min fairness. 
Table[T]shows our measurements for all instances compOl, comp02, . . . , comp21 from the ITC2007 
competition (see [IDGS12II for instance data). Please note that the known best timetables were not 
created with fairness in mind, but the objective was to create timetables with minimal total penalty. 
In Table [H Sbest refers to the known best solution for each instance. A and A' refer to the allocation 
functions given in ^ and dH), respectively. The data indicates that the timetables with a low total 
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Table 2: The performance of MaxMinFair_S A with AE = A^cw for different values of d. 



5 


10" 


10-2 


10-^ 


10-^ 


10^ 




02 


02,05 


02 


10-^ 


10,19,20 




09 


19 


10--^ 


01,10,19,20 






03,19 


10"^ 


01,10,20 









penalty are also rather fair. This can be explained by the fact that these timetables do not have a 
large amount of penalty to distribute over the curricula. Thus, most curricula receive little or no 
penalty and consequently, the distribution is fair for most curricula. We will show in Section 15.21 
however, that for timetables with a comparatively large total penalty there is still some room for 
improvement concerning fairness. 

The rightmost column of Table [T] contains the sorted allocation vectors of the best solutions. 
For a more convenient presentation, all entries of the sorted allocation vectors are multiplied by -1. 
The exponents denote how often a certain number occurs. For example, the sorted allocation vector 
(—5, —5,0,0,0) would be represented as 5^,0^. The sum of the values of an allocation vector is 
generally much larger than the total penalty shown in the second column. The reason for this is 
that the penalty assigned to a course is counted for each curriculum the course belongs to. With a 
few exceptions the general theme seems to be that the penalty is assigned to only a few curricula 
while a majority of curricula receives no penalty. In the next section we will show that the situation 
for the curricula which receive the most penalty can be improved with max-min fair optimization 
for 15 out of 21 instances. 



5.1 Max-Min Fair Optimization 

In Section m we presented algorithm MaxMinFair_SA for solving max-min fair minimization 
problems. A crucial part of this algorithm is the energy difference measure which determines how 
much worse a given solution is compared to another solution, i.e. the energy difference of the 
solutions. We evaluate the impact of the three energy difference measures (|5]), (|7]) and ^ on the 
performance of MaxMinFair_SA. 

Our test setup was the following: For each energy difference function we independently per- 
formed 50 runs with MaxMinFair_S A. The temperature levels were determined experimentally, 
we set -drnax = 5 and '6-min = 0.01; the cooling rate a was set according to ©. In order to establish 
consistent experimental conditions for fair optimization, we used a timeout, which was determined 
according to the publicly available ITC2007 benchmark executable. On our machines (i7 CPUs 
running at 3.4GHz, 8GB RAM), the timeout was set to 192 seconds. The MaxMinFair_SA 
algorithm was executed on a single core. We generated feasible initial timetables for MaxMin- 
Fair_SA as a preprocess using sequential heuristics proposed by [jBMMiQTl. The soft constraint 



violations were not considered at this stage. Since the preprocess was performed only once per 
instance (not per run), it is not counted in the timeout. However, the time it took was negligible 
compared to the timeout (less than 1 second per instance). 
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Table 3: The performance of MaxMinFair_SA with energy difference measures tsEiex, ^Eps and 





















all except 01,06,08,11,17 




18 


AEcw 


all except 1 1 


06,07,08,17,21 





Table [21 shows the impact of the parameter 5 on the performance of MaxMinFair_SA with 
energy difference measure AEcw- For each pair of values we performed the one-sided Wilcoxon 
Rank-Sum test with a significance level of 0.01. The data indicates that for best performance, 5 
should be small, but not too small. For 5 = 1, MaxMinFair_SA beats the other shown configu- 
rations on instance comp02 but performs worse than the other configurations on instances complO 
and comp20. For d = 10^ the overall performance is better than for 5 = 1, but worse than for the 
other configurations. With 5 = 10^^ and 5 = 10^^, MaxMinFair_SA shows the best relative 
performance. Thus, for our further evaluation we set 5 = 10 

Table [3] shows the relative performance of Algorithm MaxMinFair_SA for the proposed en- 
ergy difference measures (|5}, <^} and dS). The table shows for any choice of two energy difference 
measures / and j, for which instances MaxMinFair_SA with measure i performs significantly 
better than MaxMinFair_SA using measure j. Again, we used the Wilcoxon Rank-Sum test with 
a significance level of 1 percent. The data shows that AE^-w is the best choice among the three alter- 
natives, since it is a better choice than AEjex on all instances except compl 1 and a better choice than 
AEps on five out of 21 instances. However, although AEcw shows significantly better performance 
than the other energy difference measures, it did not necessarily produce the best timetables on all 
instances. For the instances compOS, compl5, comp05 and compl2 for example, the best solution 
found with AE = AEps was better than with AE = AEcw 

The data in Table |4] shows a comparison of the sorted allocation vectors of the known best 
solutions from the CB-CTT website by [DGS12J with the best solutions found by the 50 runs of 
MaxMinFair_SA with AE = AEcw First of all, for instances compOl and compl 1, the allocation 
vectors of the best existing solutions and the best solution found by MaxMinFair_SA are iden- 
tical. This means that MaxMinFair_SA finds reasonably good solutions despite the certainly 
more complex fitness landscape due to max-min fair optimization. We can also observe that the 
maximum penalty any curriculum receives is significantly less for most instances and the penalty 
is more evenly distributed across the curricula. This means that although max-min fair timetables 
may have a higher total penalty, they might be more attractive from the students' perspective, since 
in the first place each student notices an unfortunate arrangement of his/her timetable, which is tied 
to the curriculum. Furthermore, we can observe that if the total penalty of a known best solution is 
rather low, then it is also good with respect to max-min fairness. For several instances in this cate- 
gory, (compOl, comp04, comp07, complO and comp20), the solution found by MaxMinFair_SA 
is not as good as the known best solution with respect to max-min fairness. We can conclude that 
if there is not much penalty to distribute between the stakeholders, it is not necessary to enforce a 
fair distribution of penalty. 
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Table 4: Comparison of the sorted allocation vectors of the known best solutions from the website 
by [ DGS12II with the allocation vectors found by MaxMinFair_SA with respect to max-min 
fairness. 



Instance 


Known best solution 


MaxMinFair.SA (A£ = ^Ec„) 


compOl 




52,012 


comp02 




42,231,1', 03" 


comp03 


l3,lo^9,7^6^5l^4,2^o3■' 


^4 4II 222 l3 028 


comp04 


7,63, 54,42, 2,0^" 


64,4^,24, 1,046 


comp05 


4l2,36^35^32^31^30^28,...,2,03 


192 , 18^ , , 16^ , 152 , 14^5 , 13^ . . , 4^ 3^ 2 


comp06 


12,72,54,23,0*^° 


12,42,230,113,024 


compOT 


6,0'" 


6,223,124,029 


compOS 


7,63,54,42,22,0^9 


64,42,2', 1^,043 


comp09 


105,9,710, 6^5lO, 4,2,041 


6',4i4,2i',035 


complO 


22 q65 


219,1^,042 


compll 


013 


013 


compl2 


45, 3014 , 28 , 272 , 26^ , 25 1^ 224 , . . . , 22 , 03 


103,9*,83i,7',643,52,43'*,32,2i*,l,03 


complS 


8,7,6^,5^42,23,04' 


6*,44,2i3,l6,03' 


compl4 


84,7,52, 2^ 04' 


84,42,3,2l^035 


compl5 


103,93, 7, 64, 513,4,2', 03*^ 


^4 4II 223 l2 028 


compl6 


72, 5', 4,0^1 


45 2I6 l4 o46 


complT 


102,63,5^24,052 


102,62,4',3,22S,l',026 


complS 


17,15,14,13,11,10,92,51^22,023 


420 211,15 0I6 


compl9 


13,7,64, 52, 4,2', O^o 


64,4^2l^ll4,o2' 


comp20 


22 q76 


45,33,231,1', 032 


comp21 


12, 11, 104,9,74,64,512,4,23, l2,045 


10,64,5,41^3,236, 13,0^' 
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5.2 The Tradeoff Between Fairness and Efficiency 



We proposed the JFI-CB-CTT problem formulation in Section |3l which allows us to investigate 
the tradeoff between fairness and efficiency which arises in course timetabling. We can observe 
in column 4 of Table [U that for all of the best solutions from [IDGS121I the fairness index (H) 
is greater than 0.8, i.e., the known best solutions are also fair for more than 80 percent of the 
curricula. In order to solve the corresponding JFI-CB-CTT instances, we use the multi-objective 
optimization algorithm AMOSA proposed by [BSMD08J that is based on simulated annealing like 
Algorithm MaxMinFair_S A. Since we do not expect from a general multi-objective optimization 
algorithm to produce solutions as good as the best CB-CTT solvers, we will consider the following 
scenario to explore the tradeoffs between fairness and efficiency: starting from the known best 
solution we examine how much increase in total penalty we have to tolerate in order to increase 
the fairness further. We will take as examples the six instances with the highest total amount of 
penalty, compOS, comp05, comp09, compl2, compl5 and comp21. 

The temperature levels for the AMOSA algorithm were set to i^max = 20 and T^mm = 0.01; a 
was set according to @ with a timeout determined by the official ITC2007 benchmark. The plots 
in Figure[T]show the (Pareto-) non-dominated solutions found by AMOSA. The arrows point to the 
starting point, i.e. the best available solutions to the corresponding instances. For instances comp05 
and comp21 solutions with a lower total cost than the the previously known best solutions were 
discovered by this approach. The plots show that the price for increasing the fairness is generally 
not very high - up to a certain level, which depends on the instance. In fact, for comp09 and 
comp21, the fairness index can be increased by 3.5 percent and 1.4 percent, respectively, without 
increasing the total penalty at all. 

In Figure [H the straight lines that go through the initial solutions show a possible tradeoff 
between fairness and efficiency: the slopes were determined such that a 1 percent increase in 
faimess yields a 1 percent increase in penalty. For the instances shown in Figured] the solutions 
remain close to the tradeoff lines up to a faimess of 94 to 97 percent, while a further increase 
in faimess demands a significant increase in total cost. For the instances comp05, comp09 and 
compl5, there are several solutions below the tradeoff lines. Picking any of the solutions below 
these lines would result in an increased faimess without an equally large increase in the amount 
of penalty. This means picking a fairer solution might well be an attractive option in a real-world 
academic timetabling context. For compOS for example, the faimess of the formerly best known 
solution with a total penalty of 291 can be increased by 5.4 percent at 302 total penalty, which is a 
3.8 percent increase. 

In summary, improving the faimess of an efficient timetable as a post-processing step seems 
like a viable approach for practical decision making. Using a very efficient solution as a starting 
point means that we can benefit from the existing very good approaches to creating timetables 
with minimal total cost and provide improved faimess depending on the actual, instance-dependent 
tradeoff. 
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Figure 1: Non-dominated solutions found by the AMOSA algorithm for the JFI-CB-CTT versions 
of instances compOS, comp05, comp09, compl2, compl5 and comp21. All graphs show the fairness 
index on the horizontal axis and the amount of penalty on the vertical axis. 
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6 Conclusion 



In this paper we introduced two new problem formulations for academic course timetabling based 
on the CB-CTT problem model from track three of the ITC2007, MMF-CB-CTT and JFI-CB- 
CTT. Both problem formulations are aimed at creating fair course timetables in the setting of a 
university but include different notions of fairness. Fairness in our setting means that the penalty 
assigned to a timetable is distributed in a fair way among the different curricula. The MMF- 
CB-CTT formulation aims at creating max-min fair course timetables while JFI-CB-CTT is a bi- 
objective problem formulation based on Jain's fairness index. The motivation for the JFI-CB-CTT 
formulation is to explore the tradeoff between a fair penalty distribution and a low total penalty. 

Furthermore, we proposed an optimization algorithm based on simulated annealing for solving 
MMF-CB-CTT problems. A critical part of the algorithm is concerned with measuring the en- 
ergy difference between two timetables, i.e., how much worse a timetable is compared to another 
timetable with respect to max-min fairness. We evaluated the performance of the proposed algo- 
rithm for three different energy difference measures on the 21 CB-CTT benchmark instances. Our 
results show clearly that the algorithm performs best with A^cw as energy difference measure. 

Additionally, we investigated the fairness of the known best solutions of the 21 CB-CTT in- 
stances with respect to max-min fairness and Jain's fairness index. These solutions were not created 
with fairness in mind, but our results show that all of the solutions have a fairness index greater 
than 0.8. This means they can be considered quite fair. Nevertheless, our results show that some 
improvements are possible with respect to both max-min fairness and Jain's fairness index. The 
timetables produced by our proposed MaxMinFair_SA algorithms are better than the known 
best ones with respect to max-min fairness for 15 out of 21 instances. Our investigation of the 
tradeoff between fairness and the total amount of penalty using the JFI-CB-CTT problem formula- 
tion shows that the fairness of the known best timetables can be increased further with only a small 
increase of the total penalty. 
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